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Abstract 

We show that the present experimental LEP average R = I\/T/ = 20.795 ± 0.040 is 
not unambiguous due to the presence of substantial systematic effects which cannot be 
interpreted within gaussian statistics. We find by Montecarlo simulation that the C.L. 
of the original LEP sample is only 3.8 ■ 10~ 4 . We suggest that a reliable extimate of the 
true R-value is 20.60 < R < 20.98 which produces only a very poor determination 
of the strong coupling constant at the Z mass scale, 0.10 < a s (M z ) < 0.15. 
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The determination of the strong coupling constant at the Z-mass scale a s (M z ) is of pri- 
mary importance for a consistency check of perturbative QCD. In this context, the quantity R, 
defined as the ratio between the hadronic and the leptonic partial widths of the Z boson, plays a 
fundamental role. Indeed, this particular observable, operatively defined through the ratio of the 
peak cross-sections in the corresponding hadronic and leptonic channels, can determine a s (M z ) 
to a very high degree of accuracy thus allowing a direct comparison with the perturbative evo- 
lution of a s from precise low-energy data for Deep Inelastic Scattering (DIS). 

The theoretical prediction at one-loop in the electroweak theory and including O(a^) 
perturbative QCD corrections, can be conveniently expressed by using the result of the recent 
analysis by Hebbeker, Martinez, Passarino and Quast [[lj] as 

R Th = i*M (1 + 6 QCD ) (1) 

where R(°> is the purely electroweak value in the quark-parton model and 5qcd is conveniently 
expressed as [[]]] 

5 QCD = 1.06^ + 0.9(^) 2 -15(^) 3 (2) 
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By using the experimental LEP average presented at the Glasgow Conference [[|] 

R LEP = 20.795 ± 0.040 (3) 

one deduces the value [@] 

a s {M z ) =0.126 ±0.006 (4) 
or, by including all lineshape data, [|3|] 

a L s EP (M z ) = 0.127 ± 0.005 (5) 

Eqs.(4,5) should be compared with the prediction |Q] from DIS (including a fair extimate 
of the theoretical error) 

af IS {M z ) = 0.113 ± 0.005 (6) 

As pointed out by Shifman [|5|], the discrepancy between Eqs.(4,5) and Eq.(6) is disturb- 
ing, implying a rather large difference in the values of the QCD scale parameter (in the MS 
scheme and with five flavours), namely Aq CD ~ 500 MeV rather than the value Aq CD ~ 200 
MeV expected from the QCD sum rules based on the Operator Product Expansion approach. 

The presence of a possible discrepancy with the low energy extrapolations provides a 
valid motivation to reconsider critically the meaning of the experimental LEP average presented 
in Eq.(3). Indeed, as we shall explicitly show in the following, the interpretation of the exper- 
imental data is not unambiguous and the average in Eq.(3) is faced with serious problems of 
statistical consistency. 

The individual LEP measurements of R in the various fi, r and electron channels, as 
presented by ALEPH, DELPHI, L3 and OPAL at the Glasgow Conference and summarized in 
ref. [@], are reported in Table 1. 

These 12 individual measurements are not all statistically independent. However, in a 
first approximation, if one neglects the small correlation among measurements in the same 
experiment and treats all i?-values as independent, one gets precisely the same average as ob- 
tained in ref. [|2|] by using the full covariance matrix. Thus, to good approximation, one may 
be tempted to consider the 12 individual measurements in Table 1 as belonging to a normal 
population governed by gaussian statistics. 
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41.59±0.13 
20.67±0.13 
20.91±0.14 
20.69±0.12 


41.26±0.17 
20.96±0.16 
20.60±0.12 
20.64±0.16 


41.44±0.15 
20.94±0.13 
20.93±0.14 
20.70±0.17 


41.47±0.16 
20.90±0.13 
20.855±0.097 
20.91±0.13 



Table 1: R; values of the four LEP experiments. 
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Figure 1: Distribution of the 12 R t determinations 

To understand the possible presence of systematic effects, which can affect the global 
average in an uncontrolled way, we started reporting in a histogram the central values of the 12 
individual measurements. 

By inspection of fig.l one discovers the following unexpected result: near the global 
average R = 20.795 ± 0.040, where there should be a very large number of data, one finds 
a minimum of the probability since no experiment, in any individual channel, is reporting a 
central value lying in the interval 20.755 — 20.835. The various measurements, instead, can be 
divided into two sets rather sharply peaked around R ~ 20.92 and R ~ 20.66. 

In order to have a better qualitative understanding of the problem we have reported the 
12 experimental data in sequence in fig. 2 with their errors. 

As one can see, all points lie at ~ la from the central value so that the x 2 is good indeed. 
However, a good value of the x 2 does not tell much on the gaussian nature of the data. 

To obtain a quantitative description of this statement we decided to test the hypothesis of 
the common belonging of the measurements to a normal population having the observed mean 
value R = 20.795 and errors like those of each individual measurement. The variable chosen as 
a probe of non-normality is the kurtosis 

7 = /V/"2 - 3 

where fx n — J(x — x) n f(x)dx. 

We have used a random number generator to produce a large number of equivalent 
copies starting from our original population of 12 measurements reported in Table 1. For any 
generated sample of 12 measurements, with their respective errors, we compute the mean R, 
the standard deviation a and the kurtosis 7. 
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Figure 2: The determination of R t for each experiment compared to the LEP average value 
(dashed band) 

The distribution of R is shown in fig. 3 for 10000 generated configurations of 12 mea- 
surements. 



I Mean 20.79 
| RMS 0.3973E-01 
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Figure 3: R from the MonteCarlo simulation of 10000 experiments in which the errors are 
assumed to be purely statistical 



As one can see the value for R and its a are equal to those of LEP measurement (3) 
confirming the substantial statistical nature of each individual error. This provides a check of 
our approximation in neglecting the possible correlations among the errors in Table 1 . 

Fig. 4, on the other hand, shows that the probability of the initial LEP configuration in 
Table 1 is extremely small. 

The MonteCarlo runs at high statistics (10 6 trials) show that the probability to have a 
result worse than the one observed in this experiment is 3.8 • 10~ 4 . The fact that the kurtosis 
distribution does not look what is expected for a gaussian population (null mean value and a 
symmetric distribution) depends on the fact that the estimator is biassed and only asymptotically 
gets to the expected value (12 samples are not close to infinity !). This circumnstance is only 
aesthetical and does not affect the point we made. 
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Figure 4: Distribution of the kurtosis values obtained from the MonteCarlo simulation of 50000 
experiments. The arrow shows where the actual LEP results fall. 

In conclusion, our analysis indicates that the individual measurements in Table 1 can 
hardly be considered as belonging to a gaussian population since substantial systematic effects 
are needed to understand the kurtosis distribution in figure 4. As a consequence, the meaning 
of the global average in Eq.(3) (and therefore of Eqs.(4,5) ) is not entirely clear. This makes, 
at best, awkward a safe estimate of the true experimental i?-value from the data reported in 
Table 1. The large probability contents for R ~ 20.92 and R ~ 20.66 (see fig. 1 ) suggest that a 
reliable determination requires to define the error from the spread of the central values in fig. 2, 
i.e. from a full region 

20.60 < R < 20.98 
This range, by itself, allows only a very poor determination of a s 

0.10 < a s (M z ) < 0.15 

of comparable precision to that attainable from the total and hadronic Z widths and very far 
from the expected level of precision for LEP experiments. 
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